The rp-Process in Neutrino-driven Winds 
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ABSTRACT 

Recent hydrodynamic simulations of core-collapse supernovae with accurate neutrino trans- 
port suggest that the bulk of the early neutrino-heated ejecta is proton rich, in which the pro- 
duction of some interesting proton-rich nuclei is expected. As suggested in recent nucleosynthesis 
studies, the rapid proton-capture (rp) process takes place in such proton-rich environments by 
bypassing the waiting point nuclei with the /?+-lives of a few minutes via the faster capture of 
neutrons continuously supplied from the neutrino absorption by protons. In this study, the nucle- 
osynthesis calculations are performed with the wide ranges of the neutrino luminosities and the 
electron fractions {Ye), using the semi-analytic models of proto- neutron star winds. The masses 
of proto-neutron stars are taken to be 1.4 Mq and 2.0 Mq, where the latter is regarded as the test 
for somewhat high entropy winds (about a factor of two). For Yg > 0.52, the neutrino-induced 
rp-process takes place in many wind trajectories, and the p-nuclei up to A ~ 130 are synthesized 
with interesting amounts. However, ^^Mo is somewhat underproduced compared to those with 
similar mass numbers. For 0.46 < < 0.49, on the other hand, ^^Mo is significantly enhanced 
by the nuclear flows in the vicinity of the abundant ^'^Zi that originates from the a-process at 
higher temperature. The nucleosynthetic yields are averaged over the ejected masses of winds, 
and further the Ye distribution predicted by the recent hydrodynamic simulation of a core-collapse 
supernova. Comparison of the mass-yg-averaged yields to the solar compositions implies that the 
neutrino-driven winds can be potentially the origin of light p-nuclei up to A ~ 110, including 
92,94]yjQ g^j-^j 96,98j^^ -j-j^g^^ cannot be explained by other astrophysical sites. This parametric 
study suggests that the neutrino-induced rp-process takes place in all core-collapse supernovae, 
and likely in coUapsar jets or disk winds formed around a black hole. Although the quantitative 
prediction of rp-processed abundances requires future detailed multi-dimensional hydrodynamic 
simulations with accurate neutrino transport, the current investigation will serve unique con- 
straints to the fluid dynamics of the early supernova ejecta. 



Subject headings: nuclear reactions, nucleosynthesis, abundances - 
— supernovae: general 



- stars: abundances — stars: neutron 



1. Introduction 

The rapid proton-capture (rp) process of nucle- 
osynthesis is expected to take place in proton-rich 
compositions with sufficiently high temperature 
(Wallace & Woosley 1981), which leads to the pro- 
duction of proton-rich isotopes beyond iron. This 
is an analogous phenomenon to the rapid neutron- 
capture (r) process in neutron-rich compositions, 



which accounts for the production of about half 
of the species heavier than iron (see Wanajo & 
Ishimaru 2006b, for a recent review). A major 
difference between these two processes is the pres- 
ence of the Coulomb barrier for proton capture 
reactions. This limits the requisite temperature 
for this process to a small range, ~ 1 — 3 x 10^ K, 
in which the proton capture proceeds during short 
time compared to the relevant /J^-decay lifetimes 
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without substantial photodisintcgrations. To date 
the thermonuclear explosions of hydrogen-rich ma- 
terial on accreting neutron stars (X-ray bursts) 
have been considered to be the promising astro- 
physical site associated to the occurrence of re- 
process (e.g., Schatz et al. 1998; Koike et al. 1999; 
Brown et al. 2002; Woosley ct al. 2004). However, 
little amount of the nucleosynthetic product is ex- 
pected to be ejected, whose contribution to the 
Galactic chemical evolution is likely to be negligi- 
ble. 

Recent hydrodynamic studies of core-collapse 
supernovae that take the accurate neutrino trans- 
port into account have shown that the bulk of 
the neutrino-heated ejecta during the early phase 
(< 1 s) is proton-rich (Liebcndorfer et al. 2003; 
Pruet et al. 2005; Buras et al. 2006; Kitaura, 
Janka, & Hillebrandt 2006; Frohlich et al. 2006a). 
Given this proton richness in the innermost ejecta 
is the general characteristic, one might think that 
core-collapse supernovae provide suitable physi- 
cal conditions for the Ji^-process. If this is true, 
the study of r7>-process in core-collapse supernovae 
will be of special importance, which no doubt eject 
the nucleosynthetic products and thus contribute 
to the Galactic chemical evolution. However, there 
are a number of nuclei with the /3+-lives of a few 
minutes on the rp-process path. These "wait- 
ing point" nuclei inhibit the production of heavy 
proton-rich nuclei beyond the iron-group in core- 
collapse supernovae (Pruet et al. 2005). The first 
waiting point nucleus encountered soon after the 
abundant ^^Ni is ^^Ge. The half hfe of its /?+- 
decay is 1.06 min, which is obviously longer than 
the dynamic timescale of the innermost ejecta of 
a core-collapse supernova (< Is). 

The situation changes dramatically, however, 
when neutrino-induced reactions are included in 
the nucleosynthesis calculations. In fact, recent 
studies by Frohlich et al. (2006a,b), Pruet et al. 
(2006), and Wanajo (2006) have shown that the 
/S^-waiting points are bypassed via neutron cap- 
ture reactions even in proton-rich environments. 
This is due to the continuous supply of neutrons 
from the anti-electron neutrino absorption by free 
protons in the early ejecta that is subject to an 
intense neutrino flux. As a consequence, the re- 
process takes place, which leads to the produc- 
tion of proton-rich nuclei beyond the iron group. 
In particular, the production of some light p- 



nuclei by this neutrino-induced rp-proccss, or "i/p- 
process" , is anticipated in Frohlich et al. (2006b) 
and Pruet et al. (2006). In their works, however, 
only a small number of thermodynamic trajecto- 
ries were considered for the nucleosynthesis calcu- 
lations, which were also limited for the very early 
ejecta (^ 1 s after core bounce). This is due to 
the limitations of the hydrodynamic calculations 
with accurate neutrino transport currently avail- 
able, which led to "successful" supernova explo- 
sions. It is obvious, however, that more investiga- 
tions for various physical conditions are needed for 
a profound understanding of the overall picture of 
this newly discovered nucleosynthesis process. 

In this paper, therefore, the neutrino-induced 
rp-process of nucleosynthesis in neutrino-driven 
winds is investigated in some detail, for various 
physical conditions obtained from (much simpler) 
semi-analytic calculations. A special attention is 
paid for the production of light p-nuclei including 
92,94]^Q and ®^'®^Ru, which cannot be explained 
even by the most successful scenario (i.e., the 
0/Ne layers in core-collapse supernovae, Prant- 
zos et al. 1990; Rayet et al. 1995). The produc- 
tion of some light p-nuclei ('^''Se, "^^Kr, ^^Sr, and 
^^Mo) in slightly neutron-rich winds has also been 
suggested by Hoffman et al. (1996). Hence, the 
nucleosynthesis calculations are performed for the 
wide range of initial electron fractions (Yg) cover- 
ing both proton- rich and neutron-rich winds. This 
is of particular importance to discuss whether the 
neutrino- driven winds can be the major produc- 
tion site of these light p-nuclei. 

The thermodynamic and hydrodynamic trajec- 
tories of the neutrino-heated ejecta are obtained 
with the semi-analytic models of neiitrino- driven 
winds (§ 2), which were originally developed for 
the study of r-process (Wanajo et al. 2001). This 
study is thus a natural extension of the previ- 
ous nucleosynthesis study in Wanajo et al. (2001), 
to take the variation of into consideration 
(see also Wanajo et al. 2002; Wanajo & Ishimaru 
2006b). The masses of proto-neutron stars are 
taken to be 1.4 M© and 2.0 M© as in Wanajo et 
al. (2001). The latter is regarded as the test case 
for (reasonably) high entropy winds (about a fac- 
tor of two), which has been proposed as a pos- 
sible model for the production of heavy r-nuclei 
(Otsuki et al. 2000; Wanajo et al. 2001). The 
nucleosynthesis in each wind trajectory is cal- 
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culatcd for the wide range of Ye including both 
proton-rich and neutron-rich compositions (§ 3). 
The mass-averaged yields over the various neu- 
trino luminosities for a given initial arc com- 
pared with the solar abundances (§4). These re- 
sults are further averaged over the Ye distribution 
of the ncutrino-hcatcd cjccta obtained from the 
recent two-dimensional hydrodynamic calculation 
by Buras et al. (2006). These mass-Fe-averaged 
yields are compared to the solar compositions to 
discuss the contribution of the nucleosynthesis in 
neutrino- driven winds to the Galactic chemical 
evolution. Summary and conclusions of the cur- 
rent study are presented in § 5. 

2. Neutrino-driven Wind Models 

After several 100 ms from the core bounce, hot 
convective bubbles are evacuated from the surface 
of a proto- neutron star, and the winds driven by 
neutrino heating emerge, as can be seen in some 
hydrodynamic simulations of "successful" super- 
nova explosions (Woosley et al. 1994; Buras et al. 
2006). During this wind phase, a steady flow ap- 
proximation may be justified. The wind trajecto- 
ries in this paper are thus obtained using the semi- 
analytic models of neutrino-driven winds in Ot- 
suki et al. (2000) and Wanajo et al. (2001, 2002), 
which were developed for the study of r-process 
(see also Qian & Woosley 1996; Cardall & Fuller 
1997; Thompson, Burrows, & Meyer 2001). Here 
the models and some modifications added to them 
are briefly described. 

In this study, a couple of hydrodynamic studies 
of "exploding" core-collapse supernovae are closely 
referred, in order to link the current parametrized 
models to more realistic ones. One is the one- 
dimensional core-collapse simulation of a 20 Mq 
progenitor star in Woosley et al. (1994). This is in 
fact only the case that followed the wind phase for 
long duration (up to ~ 20 s), which is relevant for 
the current study. It is cautioned, however, that 
the neutrino transport was treated rather approx- 
imately, which makes it difficult to obtain the reli- 
able time evolution of Y^ in the neutrino-processed 
material. Note also that the spatial Y^ distribution 
cannot be obtained with one-dimensional simula- 
tions. The other is the two dimensional simula- 
tion of a collapsing 15 Mq progenitor by Buras 
et al. (2006), which is also taken in Pruet et al. 



(2005, 2006). The simulation was limited to the 
first ~ 0.5 s after core bounce (and linked to a one- 
dimensional simulation up to ~ 1.3 s), but the ac- 
curate neutrino transport was taken into account. 
Therefore, the obtained time variation of Yg and 
its spatial distribution may be reliable. 

The system is treated as time stationary and 
spherically symmetric, and the radius of the neu- 
tron star is assumed to be the same as that of 
the neutrino sphere. The physical variables in the 
wind are then expressed as functions of the dis- 
tance r from the center of the neutron star. The 
ejected mass by neutrino heating is assumed to 
be negligible compared to the mass of the neutron 
star. Therefore, the gravitational field in which 
the neutrino-heated matter moves can be treated 
as a fixed- background Schwarzschild metric. The 
velocity v{r), temperature T(r), and density p{r) 
can be solved with their boundary conditions by 
use of relations based on baryon, momentum, and 
energy conservation (eqs. (l)-(3) in Wanajo et al. 
2001). The source term in the equation of energy 
conservation is due to both heating and cooling by 
neutrino interactions. The gravitational redshift 
of the neutrino energies, and the bending of the 
neutrino trajectories expected from general rela- 
tivistic effects, are explicitly taken into account 
(Otsuki et al. 2000). The neutrino luminosities 
Li, of all neutrino flavors are assumed to be equal, 
and the rms average neutrino energies are taken to 
be 10, 20, and 30 MeV, for electron, anti-electron, 
and the other flavors of neutrinos, respectively. 

As the boundary conditions at the neutrino 
sphere, the density is taken to be 10^" g cm~^ and 
the temperature is determined so that heating and 
cooling by neutrino interactions arc in equilibrium. 
The mass ejection rate at the nciitrino sphere M 
that determines the initial velocity is set to be 
M = 0.99 X Mc- Here Mc is the critical value that 
gives the transonic solution. The wind with Mc 
becomes supersonic through the sonic point, while 
those with M < Mc are subsonic throughout. In 
fact, recent hydrodynamic simulations show that 
the fast wind collides with the dense shell of slower 
ejecta behind the shock and is decelerated again to 
be subsonic (e.g., Buras et al. 2006). On the other 
hand, the subsonic wind, which cannot escape by 
itself from the gravitational potential, must be 
enough fast (i.e., M w Mc) to reach the early 
bubble ejecta. For example, the wind is too slow 
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to reach the ejecta behind the shock (at several 
1000 km) when M is taken to be, e.g., 0.9 x Mc. 
Thus, the choice of M above may be reasonable. 

In this study, the neutron star masses M are 
taken to be 1.4 Mq and 2.0 Mq. The radius 
of neutrino sphere is assumed to be R^{L^) = 
{Ri,o — i?^i)(L^/L^o) + Ri^i as a function of L^, 
where R^^o — 30 km, R^i = 10 km, and L^^q 
■1 X 10^^ ergs s^^, as shown in Figure 1 {dashed 
line). This mimics the early time evolution of 
neutrino sphere that can be seen in realistic core- 
collapse simulations (Woosley et al. 1994; Buras 
et al. 2006, see also Figures 2 and 3). In Fig- 
ure 1, M for each M is also shown as a func- 
tion of Li, {solid line). The resulted (asymptotic) 
entropy per baryon s and the dynamic timescale 
Tdyn (= \p/{dp/dt)\T=o.5Mcv) are denoted by thick 
solid and dot dashed lines, respectively. As can be 
seen, the model of M = 2.0 Mq results in higher s, 
shorter Tdyn, and smaller M. In particular, the en- 
tropy is about as twice as that for the M = 1.4 Mq 
model all the way, which might lead to the pro- 
duction of the heavy r-process nuclei (Otsuki et 
al. 2000; Wanajo et al. 2001). 

It should be noted that no hydrodynamic stud- 
ies to date indicate the formation of such a massive 
proto- neutron star (with the radius of ~ 10 km). 
On the other hand, a recent hydrodynamic sim- 
ulation shows an entropy increase owing to the 
collision with the slower preceding ejecta, result- 
ing s ~ 80 at early times (~ Is, Pruet et 
al. 2005; Buras ct al. 2006). Furthermore, some 
other mechanisms that increase entropy have been 
also proposed (e.g., Qian & Woosley 1996; Thomp- 
son 2003; Suzuki & Nagataki 2005; Burrows ct al. 
2006). In fact, when it comes to merely entropy, 
the model of M = 2.0 is rather close to the 
hydrodynamic implication in Pruet et al. (2005) 
as far as during the early epoch. Thus, this model 
can be regarded simply as a reasonably high en- 
tropy case that would be expected in realistic sim- 
ulations, rather than an extreme case with a very 
massive proto-neutron star. This may be better 
than to obtain the higher entropy by multiplying 
a certain factor to the density or temperature with 
no physical basis. It should be stressed that the 
time variations of T, p, and r are obtained consis- 
tently in the current study, which are all important 
for the nucleosynthesis (§ 3.3). 

The wind trajectories are calculated for 54 con- 



stant Li, of 40 — 0.5 X 10^-^ ergs s ^ with the inter- 
vals of 1 X 10^^ ergs s~^ and 1 x 10^" ergs s~^ for 
L^>2x 10^1 ergs s'^ and < 2 x 10^^ ergs s-\ 
respectively. The set of time stationary solutions 
v{r), T{r), and p{r) are then converted to the time 
variables r{t), T{t), and p{t) for the nucleosynthe- 
sis calculations, where t is the time that is set to 
zero at the neutrino sphere for a given wind. 

In order to link the series of constant L^, tra- 
jectories to realistic timci-evolving winds, the post 
bounce time for each L^, is defined as tpb{L,,) = 
to {L^/L^f))~^ , where to = 0.2 s (Figure 1, dot- 
ted line). This is needed to determine the ejected 
mass by each wind (i.e., M{L,,) Atpb) when cal- 
culating the mass-averaged nucleosynthesis yields 
(§4). This gives L^{tph) = L^o{tph/to)^^ and 
Ru{tpb) = {Ruo - Rui){tph/to)-^ + R.1 (Figs. 2 
and 3), which approximately mimic the hydrody- 
namic results in Woosley et al. (1994). 

The time evolutions of r, Tg (= T/IO^K), 
p, and s (specific entropy) for selected (odd- 
numbered) trajectories are shown in Figures 2 (for 
M = 1.4 Mq) and 3 (for M = 2.0 Mq). In these 
figures, the time is set to i = at tp^ = tp\,{Li,) 
for each wind. As can be seen in the top panels 
of these figures, the time evolutions of radii for 
^pb > 0.5 s are very similar to those found in the 
hydrodynamic results (Woosley et al. 1994; Buras 
et al. 2006), at which the current wind solutions 
may be justified. On the other hand, the early 
time evolutions for tp\y < 0.5 s seem less realistic, 
in which the convcctive bubbles dominate in a re- 
alistic simulation (Buras et al. 2006). This may 
not cause a severe problem, since the rp-process 
takes place dominantly at tpy, > 0.5 s as discussed 
in § 3.3. Both models (M = 1.4 Afo and 2.0 Mq) 
show qualitatively similar thermodynamic histo- 
ries as can be seen in Figures 2 and 3. Note, 
however, that the entropy is about a factor of two 
higher for the M = 2.0 Mq model all the way. 

3. Nucleosynthesis in Winds 
3.1. Reaction Network 

Adopting the wind trajectories discussed in § 2 
for the physical conditions, the nucleosynthetic 
yields arc obtained by solving an extensive nu- 
clear reaction network code. The network con- 
sists of 6300 species between the proton and neu- 
tron drip lines predicted by a recent fully micro- 
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scopic mass formula (HFB-9, Goricly ct al. 2005), 
all the way from single neutrons and protons up 
to the Z = 110 isotopes. All relevant reactions, 
i.e. (n,7), (p,7), (a, 7), (a,n), {a,p), and 

their inverse are included. The experimental data 
whenever available and the theoretical predictions 
for light nuclei {Z < 10) arc taken from the REA- 
CLIB compilation (F. -K. Thielemann 1995, pri- 
vate communication) . All other reaction rates are 
taken from the Hauscr-Fcshbach rates of BRUS- 
LUB (Aikawa et al. 2005) making use of experi- 
mental masses (Audi, Wapstra, & Thibault 2003) 
whenever available or the HFB-9 mass predictions 
(Goriely et al. 2005) otherwise. The photodis- 
integration rates are deduced applying the reci- 
procity theorem with the nuclear masses consid- 
ered. The /3-decay rates are taken from the gross 
theory predictions (GT2, Tachibana, Yamada, & 
Yoshida 1990), obtained with the HFB-9 Qp pre- 
dictions (T. Tachibana 2005, private communica- 
tion). Electron capture reactions on free nucleons 
as well as on heavy nuclei are also included (Fuller, 
Fowler, & Newman 1982; Langanke & Martinez- 
Pinedo 2001). 

Rates for neutrino capture on free nucleons 

Ue+n ^ p + e~ (1) 

De+P n + e+ (2) 

and *He, 

!/e+^He ^ ^Re+p + e- (3) 

i7e-F^He ^ ^'a + n + e+ (4) 

as well as for neutrino spallation of free nucleons 
from '*He, 

v + '^Ue ^R + p + u' (5) 

ly + '^Ue ^Rc + n + iy' (6) 

are also included (Woosley et al. 1990; McLaugh- 
lin, Fuller, & Wilson 1996). Neutrino- induced re- 
actions of heavy nuclei are not included in this 
study, which may have only minor effects (Meyer, 
McLaughlin, & Fuller 1998). The rates for neu- 
trino reactions above are expressed as (Qian et al. 
1997) 



where (-E^) and (a^,) are the average energy and 
cross section of the neutrino species responsible 
for the reaction. The time evolution of Li, (that 
is taken to be constant when deriving the wind 
trajectories) is taken into account, according to 
L,j{tpb) in § 2. The distance r from the center 
of the neutron star for each wind as a function of 
time can be seen in Figures 2 and 3. 

3.2. Initial Compositions and 1^ 

Each calculation is initiated when the temper- 
ature decreases to Tg = 9 (where Tg = T/IO^K). 
At this high temperature, the compositions in the 
nuclear statistical equilibrium (NSE) are obtained 
(mostly free nucleons and a particles) immedi- 
ately after the calculation starts. The initial com- 
positions are thus given by Xn = 1 — Yei and 
Xp = Ygi, respectively, where X„ and Xp are 
the mass fractions of neutrons and protons, and 
Yei is the initial electron fraction (muiibcr of pro- 
ton per nucleon) at Tg = 9. In a core-collapse 
simulation for long duration (~ 20 s, Woosley et 
al. 1994), it is observed that the initial Yg takes 
mostly a constant value (« 0.46) during the early 
phase {tpb = a few s), and gradually decreases to 
w 0.38. This is due to the hardening of the anti- 
electron neutrino spectrum relative to that of the 
electron neutrino at later times (Qian & Woosley 
1996). In order to mimic this effect, Yei is as- 
sumed to be constant {Ygo) for to < tpb < ti 
and Fei(ipb) = (^o - i;i)(ipbAi)-^ + l^i for 
tph > ti, where ii = 4 s and Yei = 0.1. This 
gives Yei for each wind such as Y^i ~ Ygo and 
Fei(i.) = {¥,0 - Yei){L,/L,o){h/to) + Fei for 
L„>2x 10^1 ergs s'^ and < 2 x lO^^ ergs s-\ 
respectively. 

The core-collapse simulation in Woosley et al. 
(1994) shows that the neutrino-heated ejecta at 
the first few seconds is neutron rich {Yg « 0.46). 
On the other hand, recent hydrodynamic studies 
with more accurate neutrino transport show that 
the bulk of the ejecta during the early phase is 
proton rich (Y^ > 0.5, Liebendorfer et al. 2003; 
Kitaura, Janka, & Hillebrandt 2006; Buras et al. 
2006). The reason of this proton richness is that, 
when the degeneracy is lifted, the mass difference 
between the proton and neutron favors the proton- 
rich composition as far as the differences of lumi- 
nosities and mean energies between the electron 
neutrinos and the anti-electron neutrinos are not 
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significant. As a result, the electron neutrino cap- 
ture (eq. [1]) and the positron capture (inverse 
of eq. [2]) overcome their inverses, resulting in 
Y,, > 0.5 (sec Hoffman ct al. 1996; Frohlich ct al. 
2006a; Buras et al. 2006, for more detailed dis- 
cussion). In this study, the nucleosynthesis cal- 
culation for each wind is repeated for various Yeo 
(= 0.45-0.65) with the interval of 0.005 (41 cases), 
covering the slightly neutron-rich to very proton- 
rich compositions. In Figure 4 {top panel), the 
time variations of Ygi that determines the initial 
composition for each wind are shown for the winds 
of M = 2.0 Mq with YeO = 0.60 (thin solid line) 
and M = 1.4 Mq with Yeo = 0.48 {thin dotted 
line), respectively. 

Neutrino capture on free nucleons changes Y^ 
{"a effect", McLaughlin, Fuller, & Wilson 1996; 
Meyer, McLaughlin, & Fuller 1998). Figure 5 
shows Ye {thick lines) along with the mass frac- 
tions {thin lines) of neutrons X„, protons Xp, a 
particles Xa, and heavy nuclei Xh {Z > 2) as 
functions of Tg . obtained from nucleosynthesis cal- 
culations for the winds with M = 2.0 Mq, L^, = 
1 X 10^2 ergs s~\ and Y^o = 0.60 (W201060, solid 
lines) and M = 1.4 M©, = 1 x lO^^ergss'S 
and Yeo = 0.48 (W141048, dotted lines), respec- 
tively^. For the wind W201060, Yg decreases sig- 
nificantly as the temperature drops. This is due to 
the increasing a particles, which continues as far 
as neutrons are supplied from the neutrino capture 
on protons (eq. [2]). The electron fraction can be 
expressed as 

ye«Xp + ^+ Yl ZY{Z,A), (8) 

Z>2,A 

where Y{Z, A) is the abundance (number fraction 
per nucleon, i.e., the mass fraction divided by A) of 
the nucleus {Z,A). For W201060, alpha particles 
dominate with a small fraction of protons as soon 
as the temperature drops to Tg ~ 7. Hence the 
second term of equation (8) governs Vg) while the 
third term is negligible. As a result, Y^ decreases 
toward 0.5. Note that the decrease of Y^ for Tg < 3 
is due to the {n,p) reactions during the rp-process 
phase and the /3+-dccays after frcczeout (§ 3.3). 
On the other hand, for W141048, a particles and 

^The wind trajectory with M = i.j Mq, = kl x 
10^^ ergs s~^, and V^o = O.mn is labeled as Wijklmn, 
hereafter. 



also heavy nuclei dominate at Tg ~ 6. Thus, Yg 
is determined by the second and third terms in 
equation (8). The sum of these two terms are close 
to the original Y^ value (= 0.48), owing to the 
slight neutron richness of heavy nuclei. As a result. 
Ye keeps its original value. 

This can be seen in Figure 4 {bottom panel), in 
which the changes of Y,,, AYg, until the temper- 
ature drops to Tg = 3 (at which the rp-process 
begins, sec § 3.3) as function of Y^q, for the winds 
with M = 1.4 M© {dotted line) and 2.0 M© {solid 
line). The change is not significant for YeO ~ 
0.47 — 0.48 as described above. For a larger Yeo 
case (> 0.5), |/\Ye| is greater for the more massive 
M {= 2.0 Mq) model at later times (tpb = 4, i.e., 

= 2x 10^^ ergs s~^). This is due to the greater 
entropy for these winds (see Figures 1-3) that leave 
more free protons for Yg > 0.5. For the same rea- 
son, \AYf, \ is also greater for a larger Ygg case. For 
a smaller YeO case (< 0.47), AYg takes a positive 
value owing to the abundant neutrons. The elec- 
tron fraction at Tg = 3, defined as Ye/ = Yei + AYe, 
is shown as a function of fpb in Figure 4 {top panel, 
thick lines) for the winds of M = 2.0 M© with 
YeO = 0.60 {solid line) and M = 1.4 M© with 
YeO = 0.48 {dotted line). At later times (fpb > 4s), 
in which the material is assumed to be rather neu- 
tron rich. Ye increases significantly owing to the 
neutrino effect and takes k, 0.34—0.36 at fpb = 16 s 
{Lj, = 5 X 10^° ergs s~^). Note that the time evo- 
lution of Yef for M = 1.4 Mq with FgO = 0.48 
(Fig. 4, top panel, thick dotted line) is similar to 
the hydrodynamic result in Woosley et al. (1994). 

3.3. The rp-Process in Proton-rich Winds 

In the proton-rich compositions with the pres- 
ence of intense neutrino flux, the proton capture 
proceeds beyond the iron group nuclei by bypass- 
ing the /S^-waiting point nuclei, as shown in re- 
cent works (Frohlich et al. 2006b, a; Pruet et al. 
2005; Wanajo 2006). This can be clearly seen 
in Figure 6, which compares the snapshots of nu- 
cleosynthesis at t = 0.22 s^ (when the tempera- 
ture drops to Tg w 2) for W141060 without {left 
panel) and with {right panel) neutrino reactions. 
Without neutrino reactions, the nuclear flow can- 
not advance beyond ^^Ge, whose /3+ half life is 

^In this section, the time is set to t = at Tg = 9. 
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1.06 minutes^. In contrast, the neutrino capture 
on abundant protons (eq. [2]) provides a fraction of 
neutrons that immediately suffer {n,p) reactions. 
Note that the {n,p) reaction acts as the same as 
the /3+-decay of a given species. As a result, the 
nuclear flow proceeds to the heavier proton-rich 
species (Fig. 6, right panel). 

It is interesting to note that the proton-to-seed 
abundance ratio Yp/Yh at this point (Figure 6, 
right panel) is still high {— 49) despite the moder- 
ate entropy (55 NAk). This is due to the slowness 
of the a + a + a — »^^C reaction that is the bot- 
tleneck for the production of heavy nuclei. This 
is in contrast to the neutron-rich compositions, in 
which the faster a + a + n ^^Be reaction fol- 
lowed by a capture is efficient. This is why the 
entropy of a few lOOiV^fc is needed for the r- 
process in neutrino- driven winds (Woosley et al. 
1994; Wanajo et al. 2001). In the proton-rich com- 
positions, however, only a moderate entropy is suf- 
ficient for the occurrence of rp-process. 

This neutrino-induced rp-process can be seen 
more clearly for the wind models with M = 
2.0 Mq. This is a consequence of their higher en- 
tropies (see Figures 1-3), which leave more pro- 
tons needed for the rp-process. Figure 7 shows 
the snapshots of nucleosynthesis calculations for 
the winds W201060 {top panels), W200865 {mid- 
dle panels), and W200265 {bottom panels) during 
(Tg « 2, left panels) and after {t w 1.0 s, right 
panels) the rp-process phase. For the same L^, 
and Yeo to W141060 (Fig. 6, right panel) but with 
M = 2.0 Mq (W201060), Yp/Yh is three times 
greater (= 151) at Tg = 2 (Fig. 7, top left), owing 
to its higher entropy (= 92 NAk). As a result, the 
nuclear flow proceeds to the Z = hO magic nuclei, 
which leads to the production of the p-nuclei with 
A ~ 110. 

The detailed nuclear flow patterns for W201060 
when the rp-process is active (Tg = 2.5) are shown 
in Figure 8. Here, the nuclear flow for each reac- 

*In the current reaction network, ®''As is omitted, whose 
proton separation energy is predicted to be negative in the 
HFB-9 mass formula (Goriely et al. 2005). In principle, 
the nuclear flow can proceed through ®®As by two proton 
capture reaction even with the (small) negative proton sep- 
aration energy, owing to the presence of the Coulomb bar- 
rier. Inclusion of this isotope would not change the current 
result, however, since the time scale of this two p reaction 
is no less than ~ 1 s (Brown et al. 2002). 



tion i — > J is deflned as 

F^^=Y{i^3)-Y{3^i) s-i. (9) 

In Figure 8, Fij and Yi (for those greater than 

10~^) are shown in logarithmic scale by arrows 
and circles, respectively. As can be seen, {p, n) 
reactions play a dominant role to carry away the 
nuclear abundances from waiting points, while the 
contribution of /3"'"-decays (off-centered arrows) is 
minor. 

The time variations of the mass fractions {thin 
lines) of neutrons X„, protons Xp, a particles Xa, 
and heavy nuclei Xh for W201060 are shown in 
Figure 9, along with those of the average mass 
number of heavy nuclei Ah and the temperature 
{thick lines), where 

Xh= ^(^'^) (10) 

Z>2,A 

and 

^''^^ E AX{Z,A). (11) 

The results with all neutrino-induced reactions 
(eqs. [l]-[6]), neutrino capture on free nucleons 
only (eqs. [1] and [2]), and no neutrino- induced 
reactions are denoted by the solid, dashed, and 
dotted lines, respectively. The initiation of the 
rp-process can be clearly seen as the sudden in- 
crease of Aft, at Tg f» 3, for the cases with neu- 
trino reactions {thick solid and thick dashed lines). 
This is due to the presence of neutrons (X„ ^ 
10~^^) by the neutrino capture on abundant pro- 
tons (eq. [2]), as can be seen in Figure 9 {thin solid 
and thin dashed lines). 

The neutrino reactions for ^He tend to inhibit 
the rp-process as can be seen in Figure 9, although 
their effects are moderate. This is due to the spal- 
lation reactions of ''He (eqs. [5] and [6]), while the 
capture reactions (eqs. [3] and [4]) are not impor- 
tant (because of their smaller cross sections). The 
resulting ^H emitted from the neutrino spallation 
of equation (5) immediately suffers a (p, n) reac- 
tion to ^He in the proton-rich compositions. These 
^He nuclei, in addition to those from the spallation 
reaction of equation (6), capture abundant a par- 
ticles to be ''Be, which further capture a particles. 
As a result, a fraction of '*He assembles into heav- 
ier nuclei that act as "neutron poison" , resulting in 
the smaller neutron abundance (Figure 9). This is 
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a similar mechanism shown by Meyer (1995), i.e., 
the r-process is also inhibited by these neutrino re- 
actions on *He in the neutron-rich compositions. 
It should be noted that the relatively large mean 
energies of /z and r neutrinos assumed in this study 
(according to Woosley et al. (1994)) may lead to 
the overestimation of this effect. 

In order to see the gross feature of the neutrino- 
induced rp-process more clearly, the time varia- 
tions of the mean lifetimes for /J^-decay, (n, 7), 
(p, 7), and {n,p) reactions, and neutrino capture 
on protons (eq. [2]) are shown in Figure 10 {solid 
lines), where 



X0+{Z,A)Y{Z,A) 



Z>2,A 



^ ^ Xr,,{Z,A)YiZ,A) 



Z>2,A 



Tpj — 



n -1 



^ ^ X^{Z,A)Y{Z,A) 

Z>2,A 

^ Xnp{Z,A)Y{Z,A) 

Z>2,A 

^ YiZ,A) 

Z>2,A 



pYn 

Yh 



and 
Here, 



(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



and Xf3+{Z,A), A„-,(Z,^), Ap^(Z,A), A„p(Z,A), 
and Ap^p are the rates of the corresponding reac- 
tions. The rates in equations (12)-(15) are aver- 
aged over the heavy nuclei, which represent the 
lifetimes of the dominant isotopes at a given time. 
The dotted lines are the corresponding lifetimes 
when the sum in equations (12)-(15) are taken for 
Z > 4:8 only. With the timcscalcs defined above, 
the requisite condition for the neutrino-induced 
rp-process can be expressed as 



(18) 



Note that the conditions for the r-process and 



classical rp-process are simply r, 
Tp^ <C Tf}+ , respectively. 



n7 



< Tfi 



and 



Top panel of Figure 10 is the result for the wind 
trajectory of W201060 that is taken to be the ref- 
erence case here. As can be seen, the {n,p) re- 
actions are faster than the /3+-decays all the way 
during the first one second. This shows clearly the 
importance of the {n,p) reactions for the neutrino- 
induced rp-process. The role of (n, 7) reactions is 
minor, compared to {n,p) reactions. A saturation 
of Ah at t = a fewlOOms (Tg ~ 2) in Figure 9 
indicates the termination of the neutrino-induced 
rp-process. At later times. Ah slightly decreases 
by photodisintegration. Note that this termina- 
tion is not due to the exhaustion of protons, as 
in the case of r-process that ceases generally by 
the neutron exhaustion (Wanajo et al. 2004). The 
proton-to-seed abundance ratio is still high at later 
times (e.g., Yp/Yh = 77 at t = Is), and the 
bulk of (p, 7) reactions are enough fast through- 
out (Fig. 10, solid line). However, the (p, 7) re- 
actions for the heaviest nuclei {Z > 48; Fig. 10, 
dotted line) become slow owing to the increasing 
Coulomb barrier, and compete with the (n,p) re- 
actions (i.e., Tpj ~ Tnp) at t ~ 0.1s (Tg ~ 2). 
This results in the broadening and smoothing of 
the abundances for Z > 40 in Figure 7 {top pan- 
els), which is similar to the "freezeout effect" seen 
in the r-process (Wanajo et al. 2004). 

What terminates the rp-process is thus the 
Coulomb barrier of the heaviest nuclei in the nu- 
clear flow, which breaks the relation between the 
flrst two in equation (18). The slowness of the 
(n,p) reactions that compete with the /3+-decays 
(t„p ~ T(3+) can be another reason to cease the 
rp-process. The (n,p) reactions becomes signifi- 
cantly slow owing to the less active neutrino cap- 
ture on protons (Figure 10) at later times {t > Is), 
whose lifetime is proportional to the square of r 
(eqs. [7] and [16]). For the current wind trajectory 
(W201060), the relation between the latter two in 
equation (18) breaks only at late times (t « 1 s for 

> 2 and t « 7s for Z > 48). This shows that 
the supply of neutrons by the neutrino capture 
on protons is sufiicient throughout the rp-process 
phase. 

One may consider that the neutrino-induced 
rp-process proceeds to much higher atomic num- 
ber to produce all the p-nuclei up to A ~ 200, 
if the proton concentration is much higher owing 
to, e.g., a larger Yeo or a higher entropy. This 
may not be true, however, as can be seen in Fig- 
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urc 7. For the wind trajectory with a larger en- 
tropy (= 103 NAk) and Yeo (= 0.65, W200865) 
that results in Yef = 0.582, Yp/Yh = 427 at Tg = 2 
and the nuelear flow proceeds beyond Z = 50 (Fig- 
ure 7, middle panels). In fact, this wind results in 
the production of the heaviest p-nuclei (A ~ 130) 
in the current study. However, the flow devi- 
ates from the proton-drip line and reaches the /3- 
stability line (Figure 7, middle right, see also Pruet 
et al. 2006). This is not due to the /3+-decays of 
the heaviest nuclei, which are rather slow all the 
way (r^+ ~ 100 s for Z > 48; Fig. 10, middle 
panel, dotted line). This is a consequence of the 
slowness of proton capture owing to the higher 
and higher Coulomb barrier for the heaviest nu- 
clei in the flow, which competes with both {n,p) 
and (n, 7) reactions as can be seen in Figure 10 
{middle panel, dotted lines). It should be also 
noted that there are a unbound nuclei loe-iosr^g 
{Z = 52), which is regarded as the end point of the 
classical ?7>-process (Schatz et al. 2001). In addi- 
tion, the proton capture also slows substantially 
when the temperature drops below Tg ~ 2 (Fig- 
ure 10). Thus, the temperature range in which the 
neutrino-induced rp-process takes place is strictly 
limited to 2 < Tg < 3. This differs significantly 
from the neutron capture without Coulomb bar- 
rier in case of the r-process. 

At the later wind phase, the entropy increases 
significantly, as can be seen in Figures 1-3. The 
situation is even worse (but interesting), however, 
for the production of proton-rich nuclei. The snap- 
shots of the nucleosynthesis calculations for the 
wind trajectory of W200265 is shown in Figure 7 
{bottom panels). This wind obtains the highest 
entropy (= 176 NAk) during the proton-rich phase 
assumed in this study (tpb < 4 s. see Fig. 4, top 
panel). As can be seen, the nuclear flow cannot 
sustain its proton richness and runs toward the /3- 
stability when the temperature drops to Tg « 2 
(Figure 7, bottom left) despite its extremely high 
yp/yh (= 858). At later times, the (n, 7) reactions 
push the flow to the neutron-rich side (Figure 7, 
bottom right). The reason is that the tempera- 
ture falls off to Tg « 1 quickly in the later winds 
(Figure 3). In addition, high Yp/Yh also results 
in (moderately) high Yn/Yh with the presence of 
intense neutrino flux. As shown in Figure 10 {bot- 
tom panel), when the temperature is below Tg ^ 2, 
the (p, 7) reactions for the heaviest nuclei in the 



flow become significantly slow, while the {n, 7) re- 
actions are still active. As can be seen in these ex- 
amples, the neutrino- induced r;p-process may not 
lead to the significant production of p-nuclei be- 
yond A ~ 130. It is interesting to note that many 
of the nuclei synthesized here (Fig. 7, bottom right 
panel) are usually assigned as s-nuclei or even r- 
nuclei, which are now synthesized in proton-rich 
compositions (see also Pruet et al. 2006). 

3.4. ^^Mo Production in Slightly Neutron- 
rich Winds 

In slightly neutron-rich compositions {Ye » 
0.47 — 0.49), two dominant quasi-statistical- 
equilibrium (QSE) clusters are formed around 
A K 60 and A « 90 during the a-process phase 
(Tg w 7 — 4). At the end of the a-process phase, 
the matter consists of mostly a particles, the 
heavy nuclei belong to these two QSE clusters, 
and a fraction of free nucleons (Figure 12). As 
soon as the a-process ceases, the proton capture 
plays an important role for the production of the 
p-isotope ^^Mo as suggested by Hoffman ct al. 
(1996), whose origin has been a long-standing 
mystery (see Arnould & Goriely 2003, for a recent 
review). The nucleosynthesis calculations covering 
the slightly neutron-rich winds are thus needed to 
estimate their possible contribution to the Galac- 
tic chemical evolution of the light p-nuclei. 

As shown in Figure 11, a fraction of the most 
dominant species ^°Zr {N = 50,Z = 40) in the 
QSE cluster (at A - 90) flows into '^^Mo {N = 
50, Z = 42) through various nuclear reaction chan- 
nels. The abundance of ^^Mo reaches 17% of its 
final value at this time (Tg = 3.5). Here, the {p, n) 
reaction plays a dominant role to carry away the 
nuclear abundances from ^^Zi. In Figure 12, the 
time variations of the mass fractions {thin lines) 
of neutrons X„, protons Xp, a particles Xa, heavy 
nuclei X^,, ^^Zr, and ^^Mo for the wind trajectory 
of W141048 are shown, along with those of the av- 
erage mass number of heavy nuclei {Ah, eq. [11]) 
and the temperature {thick lines). Note that the 
neutrino c;fFect on Yg is negligible for this wind as 
described in § 3.1, i.e., Ygf « Yeo = 0.48. The de- 
crease of Xp slows as the QSE cluster that contains 
^"Zr grows (at Tg « 5), in order to balance with 
the electron fraction of the cluster (« 0.44) that is 
smaller than that of the wind material (w 0.48). 
These protons are gradually absorbed by the nu- 
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clci in this cluster after the a-process phase. As a 
result, the abundance of ^^Mo increases, reaching 
its half final value at Tg = 3.2. 

The enhancement of ^^Mo in this way is highly 
dependent on the temperature history of the wind 
after the a-process phase. The production of 
^^Mo relics largely upon the slower (p, n) reac- 
tions (mean lifetime of about 30 ms at Tg = 3.5) 
on ^^Zr, since the photodisintcgration impedes the 
much faster (p, 7) reaction at this time (Figure 11). 
Another channel, that is the (n, 7) reaction on 
^°Zr, also contributes but suffers from its inverse. 
The (p, 7) reaction on ^"Zr becomes effective only 
when the temperature drops below Tg w 2.7 for 
W141048. Therefore, a relatively long cooling 
timescale of the wind material is required to ob- 
tain a large abundance of ^^Mo (that is close to its 
QSE value, see Hoffman et al. 1996). In the cur- 
rent study, the subsonic solutions of winds with 
M « Mf. are adopted as discussed in § 2. As a 
result, the wind is decelerated when passing over 
the sonic radius, and thus the decrease of temper- 
ature slows down. For the early winds {tph < 1 s), 
this transition takes place around Tg w 4 — 3. As 
a result, ^°Zr suffers {p,n) and (n, 7) reactions 
for longer duration, resulting in the more produc- 
tion of ^^Mo. Other light p-nuclei, e.g., ^^Mo and 
96,98j^y^ are not significantly produced in this pro- 
cess, although '''*Se and ^*Kr are moderately en- 
hanced (sec § 4, see also Hoffman et al. 1996). 

It should be noted that the flow patterns ap- 
peared in Figure 11 would be somewhat different 
when another nuclear data set were adopted (see 
Hoffman et al. 1996). In fact, most of the rele- 
vant reactions here are taken from the theoretical 
(Hauscr-Feshbach) predictions, although these nu- 
clear masses arc well determined by experiments 
(Audi, Wapstra, & Thibault 2003). Future deter- 
minations of these rates based on experiments are 
highly desirable. 

Neutrino-induced reactions arc not of impor- 
tance for the production of ^^Mo in these neutron- 
rich winds. The results with all neutrino-induced 
reactions (eqs. [l]-[6]), free nucleons only (eqs. [1] 
and [2]), and no neutrino- induced reactions are de- 
noted by the solid, dashed, and dotted lines, re- 
spectively. The dashed and dotted lines cannot be 
distinguished, which means that the neutrino cap- 
ture on free nucleons plays no role. This is due to 
their small abundances as can be seen in Figure 12. 



The neutrino capture on '*He docs not play a sig- 
nificant role, either. The spallation of free nucle- 
ons from ^He, on the other hand, slightly enhances 
the abundance of ^^Mo. The ■^H and ^He nuclei 
emitted from the spallation reactions (eqs. [5] and 
[6]) immediately capture abundant a particles to 
be ^Be and ^Li, which further capture a particles. 
The spalled free nucleons can survive owing to the 
faster a capture, resulting in the slight increase of 
^^Mo as can be seen in Figure 12. 

3.5. The r-Process in Neutron-rich Winds 

At later times {tph > 4 s), the entropy of the 
winds increases more than 100 N^k (Figs. 1-3), 
and Yei is assumed to be neutron rich (Fig. 4, top 
panel). This may provide the suitable condition 
for the production of r-process nuclei (Wooslcy et 
al. 1994; Wanajo et al. 2001). As discussed in § 4, 
the winds with M = 1.4 M0 can be the origin of 
only the light r-proccss nuclei (up to A ^ 130; 
Fig. 15, top panel). On the other hand, the heavy 
r-process nuclei (up to A ^ 190) can be produced 
in the winds with M = 2.0 Mq during the late 
phase (Fig. 15, bottom panel). This is a conse- 
quence of the higher entropies (Figs. 1-3) and in 
part the shorter dynamic timcscalcs (Fig. 1) of 
the winds, resulting in the higher neutron-to-seed 
abundance ratios at the beginning of r-processing 
(Tg ^ 3, see Wanajo et al. 2001; Wanajo & Ishi- 
maru 2006b, for more detail). Note that the nu- 
cleosynthetic results in the later winds do not con- 
tribute to the production of proton-rich isotopes 
(§4), since Yei is assumed to quickly drop to be 
rather neutron rich (Figure 4, top panel). 

4. Contribution to the Galactic Chemical 
Evolution 

4.1. Mass-averaged Yields as Functions of 

Ye 

In order to estimate the contribution of the 
neutrino-induced rp-process in neutrino-driven 
winds to the Galactic chemical evolution, the nu- 
cleosynthetic yields for each Yei model (41 cases) 
are mass-averaged over the 54 wind trajectories 
weighted by M(L,y )Zitpb. Figure 13 compares the 
averaged mass fractions Xcj of p-nuclei with re- 
spect to their solar values (Anders & Grevesse 
1989) as functions of Yg/ for M = 1.4 Mq {top 
left) and 2.0 Mq {top right) models. Here, Y^f, 
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which is approximately the value at the begin- 
ning of the rp- process phase (at Tg = 3, § 3.3), is 
taken at tph = 4s (at L,^ = 2 x 10^^ ergs s~-^) as 
representative of different Yef. The abundances 
of ^"'Zn and ^°Zr are also plotted for comparison 
purposes. The contribution of p-nuclei produced 
during the later phase (t > 4 s) is not important, 
as discussed in § 4.2. Figure 13 also shows the 
results for selected winds (not mass-averaged) for 
= 8 X 10^^ ergs s~^ (tpb = Is, middle panels) 
and I/y = 2 X 10^^ ergs s~^ {tpb = 4 s, bottom pan- 
els). The histogram that is taken from Buras et 
al. (2006) is the asymptotic Yg distribution p{Ye) 
of the neutrino-processed ejecta during the first 
468 ms after core bounce for a ISM© progenitor 
star, obtained by their two-dimensional hydrody- 
namic simulation. 

As can be seen, a variety of p-nuclei are pro- 
duced with interesting amounts for Y^j > 0.5 
models. The heavier p-nuclei (up to ^ ~ 130) 
appear for the greater Vg/ models as well as for 
the larger M (= 2.0 Mq) case (i.e., greater en- 
tropies, see Figs. 1-3). If we assume that the Y^ 
distribution by Buras et al. (2006) holds for the 
current models, the neutrino-driven winds may 
contribute the Galactic production of p-nuclei up 
to A ~ 110 C^Se, 78Kr, 84Sr, ^^'^^Mo, 96,98Ru, 
^°^Pd, and ^"^'^^^cj), which show considerable 
enhancements for Ygf « 0.50 — 0.56. For the winds 
with = 8 X 10^^ ergs s"^ (Fi gure 13, middle 
panels), it can be seen also that the heavier p- 
nuclei are highly enhanced as Y^/ increases. For 

= 2 X 10^^ ergs s~^ (Figure 13, bottom panels), 
however, each p-process abundance reaches the 
maximum and decreases again as Vg/ increases. 
The reason is that the late winds in the current 
study cool down quickly below Tg ^ 2 as can be 
seen in Figures 2 and 3. As a result, most of the 
proton-rich abundances produced during the re- 
process phase (Tg 3—2, Fig. 7, bottom left) suffer 
neutron capture to neutron-rich isotopes (Fig. 7, 
bottom right), regardless of the large Fg/ as well 
as the high entropy. 

For Yef « 0.46-0.49, some light p-nuclei (^''Se, 
^^Kr, 8^Sr, and ^^Mo) are enhanced with interest- 
ing amounts (Fig. 13). Of particular importance 
is the production of ^^Mo as discussed in § 3.4. 
Hoffman et al. (1996) showed that a neutrino- 
driven wind having Y^ > 0.484 cures the TV = 50 
overproduction (mostly ^°Zr), which is replaced 



with the production of ^^Mo. This strictly lim- 
ited Ye range is due to the competition between 
the production of ^°Zr (as the seed of ^^Mo) in 
more neutron-rich winds and the concentration 
of protons (needed for the ^^Mo production) in 
less neutron-rich winds. In the current study, the 
mass- averaged abundance of ^^Mo with respect to 
its solar value overcomes that of ^"Zr for wider 
range of Vg/ (~ 0.46 — 0.49, Fig. 13, top panels). 
The reason is that the abundances are m,ass aver- 
aged over wide range of Lj/, not from a single wind 
trajectory. In particular, the earlier winds expe- 
rience longer durations at Tg w 4 — 3, at which 
the ^^Mo abundance is greatly enhanced (§ 3.4). 
In fact, for a specific wind, the enhancement of 
^^Mo is limited to a smaller Yg/ range as can be 
seen in Figure 13 {middle and bottom panels). In 
addition, the position of Vg/ at which the ^^Mo 
abundance takes the maximum shifts toward 0.5 
for lower L^, (i.e., later winds). This is a con- 
sequence that the later winds cool down quickly 
below Tg « 4 — 3, thus the proton richness (i.e., 
larger Yg) is more important. In some later winds, 
relatively high entropies drive the material to form 
the QSE cluster that contains ''"Zr even with small 
neutron excess (Yg/ 0.49 — 0.50), resulting in the 
enhancement of ^^Mo (Fig. 13, middle right and 
bottom panels). 

The mass fractions of p-nuclei with respect 
to their solar values (Anders & Grevesse 1989) 
are also shown as functions of in Figure 14, 
for selected winds {M = 1.4 M©; left panels, 
M = 2.0 Mq] right panels, Fgo = 0.48; top pan- 
els, Ypo ~ 0.60; bottom, panels). For the winds 
with Y'eo = 0.48 (Yg/ « 0.48, as representative 
of slightly neutron-rich winds), ®^Mo is greatly 
enhanced at relatively early times {L,y > 1 x 
lO^^ergss"^, i.e., tp^ < 0.8s). This is due to 
the slow decrease of the temperature in the early 
winds at Tg « 4 — 3, in which ^^Mo is most en- 
hanced (§ 3.4). It should be cautioned that the 
steady wind approximation assumed in this study 
may not hold during the very early (bubble) phase 
{tpb < 0.5 s). However, it is likely that the ma- 
terial ejected as bubbles expand more slowly (e.g., 
Buras et al. 2006), and experiences longer time 
at Tg f» 4 — 3 in which ^^Mo might enhance to 
a similar (or more) degree. For the winds with 
Ygo = 0.60 (Yef ~ 0.54 — 0.56, as representative of 
proton-rich winds), p- nuclei are enhanced at rel- 
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atively later times (i^ < 1 x 10^^ ergs s^^, i.e., 
tph > 0.8 s), where the steady wind approxima- 
tion may hold (Figure 14, bottom panels). This is 
also consistent to the results by Pruet et al. (2006), 
which showed no significant production of p-nuclei 
in the bubbles. Note that no p-nuclei are produced 
in the very late winds (L^ < 2 x 10^"'^ ergs s"'^, i.e., 
tpb > 4 s), in which the matter is assumed to be 
rather neutron rich (Fig. 4, top panel). 

4.2. Mass-Yg-averaged Yields 

In reality, the neutrino-heated matter must 
have a certain distribution of when the multi- 
dimensional effects arc taken into account. The 
problem is, of course, the unknown mechanism of 
core-collapse supernovae, which governs the hy- 
drodynamic and thermodynamic histories of the 
neutrino-driven winds. Among a number of arti- 
ficially induced explosion models, only the two- 
dimensional calculation with accurate neutrino 
transport by Buras et al. (2006) provides us a reli- 
able Ye distribution of the neutrino-processed ma- 
terial. The explosion in their simulation was ob- 
tained by omitting the velocity-dependent terms 
from the neutrino momentum equation. This led 
to the increase of the neutrino energy deposition in 
the heating region by a few 10%, which converted a 
failed model into an exploding one. The Vg distri- 
bution p{Ye) of the neutrino-processed ejecta dur- 
ing the first 468 ms after core bounce for a 15 Mq 
progenitor star obtained by Buras et al. (2006) is 
overlaid in Figure 13. The histogram has the max- 
imum at Yg 0.5 and dominates in the proton- 
rich side. 

To test the contributions of the winds for M = 
1.4 Mq and 2.0 models, the mass-averaged 
yields for each VgO model (Fig. 13, top panels) 
arc further YJ.-avcraged 2000 winds in total for 
each M) with p{Ye), assuming this distribution to 
be representative of neutrino- driven winds. Need- 
less to say, there is no guarantee that this dis- 
tribution holds for the current wind models (at 
tpb = 4 s) for both M = 1.4 Mq and 2.0 Mq cases. 
However, this is only the case obtained by a multi- 
dimensional hydrodynamic calculation with accu- 
rate neutrino transport that is essential to predict 
the Ye distribution of the neutrino-heated ejecta. 
It may be, therefore, interesting to see the possible 
contributions to the Galactic chemical evolution 
with this p{Ye), keeping in mind that a more con- 



sistent estimation of p{Ye) will be needed in the 
future study. 

The resulting abundances with respect to their 
solar values are shown in Figure 15 for M = 
1.4 M© {top panel) and 2.0 {bottom panel) 
models as functions of mass number. The abun- 
dances smaller than X^-^/Xq < 100 are omitted 
here. The even-Z and odd-Z isotopes for a given 
element are denoted by circles and triangles, re- 
spectively, which are connected by lines. The p- 
nuclei are marked with filled symbols. The dotted 
horizontal lines indicate a "normalization band" 
(Woosley et al. 1994) between the largest produc- 
tion factor (i°"Mo for M = I.4M0 and ^^Mo for 
M = 2.0 Mq) and that by a factor of ten less than 
that, along with the median value {dashed line). 
This band is taken to be representative of the un- 
certainty in the nuclear data involved. 

When compared with the results in which the 

mass average is only for > 2 x 10^^ ergs s"'^ 
(Fig. 17, top panels), it is clear that the later winds 
for < 2 X 10^^ ergs s~^ contribute to only the 
r-nuclei with A > 90 but to p-nuclei. Figure 15 
implies that there exists a correlation between the 
production of p-nuclei and r-nuclei. That is, the 
higher entropy model (M = 2.0 M0) produces 
heavier p-nuclei (up A ~ 110) and r-nuclei (up 
to A ~ 190) in a single event. This might be 
true, since both the rp-process and r-process fa- 
vor the high entropy conditions. It would be pre- 
mature to conclude that, however, when consid- 
ering the highly uncertain late-time evolution of 
neutrino-driven winds. It is interesting to note 
that the overproduction of = 50 nuclei ^^Sr, 
^^Y, and ^°Zr seen in the previous r-process stud- 
ies (Woosley et al. 1994; Wanajo et al. 2001) arc 
now replaced with the moderate production of 
^^Mo. Figure 16 compares the results when M 
is replaced to 1.6 A^q {top panel) or 1.8 Mq {top 
panel). The higher entropy (~ 20% for the former 
and ~ 50% for the latter) with respect to the case 
of M = 1.4Mo (Fig. 15, top panel) results in the 
production of p-nuclei with higher mass numbers 
(and of r-nuclei). This shows that the moderate 
entropy of ~ 60 — 80 A^fc in early winds (~ 1 s) 
is sufficient for the reasonable production of light 
p-nuclei up to A ~ 100 - 110. 

Hereafter, only the winds for L^, > 2 x 
10^^ ergs s~^ {tpb < 4 s) are considered, in which 
the p-nuclei are produced. Note that the pro- 
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duction of Zn that is the dominant stable zinc 
isotope (but its origin is unknown, see Hoffman 
et al. 1996; Umeda & Nomoto 2002; Pruet et al. 
2005; Frohhch et al. 2006a) is about 10 times 
smaller than the lower normalization band, much 
smaller than expected in Hoffman et al. (1996, see 
also Fig. 13). The contributions to the production 
of ^^Sc and ^^Ti are not important either, which 
are suggested to be produced in the earlier ejecta 
(Pruet et al. 2005; Prohlich et al. 2006a). 

Figure 17 compares the results with {top pan- 
els) and without {bottom panels) neutrino reac- 
tions (eqs. [l]-[6]) for the M = IAMq {left pan- 
els) and 2.0 M© {right panels) models. As can 
be seen, the production of almost all ^^Mo and 
some portions of ^^Sc, ^^Kr, and *'*Sr are not due 
to the effect of neutrino reactions. This can be 
also seen in Figure 18, which compares the results 
from only neutron-rich {top panels) or proton- 
rich {bottom panels) winds for the M = IAMq 
{left panels) and 2.0 M© {right panels) models. 
^^Mo is highly overproduced in the mass aver- 
age of neutron-rich winds, which is, on the con- 
trary, somewhat underproduced in the proton-rich 
winds compared to other light p-nuclei (e.g., ^^Se, 
^^Kr, 84Sr, 94Mo, and ^^.ssru). Hence, the bulk 
of ^^Mo originates from the slightly neutron-rich 
winds {Yef w 0.46 — 0.50), while other p-nuclei 
from mainly proton-rich winds. 

When we look at the top panels of Figure 17 
again, the p-nuclei up to ^^Mo and -"^^^Cd for 
M = IAMq and 2.0 M© models, respectively, fall 
within the normalization band, which are regarded 
to be the dominant species produced by each 
event. It is interesting to note that the contribu- 
tions to °4Mo, ^f''»«Ru, and "'^pd for M = 1.4 Mq 
and "2,ii4.ii5gj.^^ ii3jj^ fyj. ^ 2.0 Mq are 

marginal. The ejected masses by winds diiring the 
first 20 s are 3.0 x IQ-^ Mq and 1.2 x 10"'^ Mq 
for the current M = 1.4 M© and 2.0 M© models, 
respectively. Given that the progenitor mass for 
each case to be, e.g., 15 M© and 30 Af©, respec- 
tively, the overproduction factor is expressed as 
~ 10~^ (Xej/X©). The requisite overproduction 
factor for the nucleosynthetic event to be the ma- 
jor source in the solar system is ~ 10 (Woosley et 
al. 1994), assuming that all the core-collapse su- 
pernovae produce the same amount of the isotope. 
This number increases (e.g., ^ 100) if a limited 
fraction of supernovae (e.g., 25— 30M©) contribute 



to these light p-nuclei production (see Ishimaru & 
Wanajo 1999; Wanajo & Ishimaru 2006b, for the 
r-process elements). This is conceivable since the 
current nucleosynthesis results are highly depen- 
dent on the physical conditions such as entropy 
or Yg. The overproduction factors of ~ 10 — 100 
(Fig. 17, top panels or Fig. 15) for the current 
models imply that the neutrino-driven winds can 
be potentially the major astrophysical site of these 
light p-nuclei. 

5. Summary and Conclusions 

In this study, the neutrino-induced rp-process 
in neutrino-driven winds has been investigated. 
The thermodynamic histories of winds were ob- 
tained from the semi-analytic models of the 
neutrino-driven winds that had been developed 
for the r-process study in previous works. The 
subsonic wind solutions with the nearly max- 
imum mass ejection rates were taken in this 
study, for various neutrino luminosities (40 — 

0. 5 X 10^^ ergs s^^) with the proto-neutron star 
masses of 1.4 M© and 2.0 Mq. The latter (2.0 Mq) 
was regarded as the test case of higher entropy 
winds (about a factor of two), which might be ex- 
pected in more realistic hydrodynamic simulations 
of core-collapse supernovae. The nucleosynthesis 
calculations were performed for the wide range of 
the initial electron fractions (0.45 < Yeo < 0.65) 
including rather proton-rich compositions, moti- 
vated by recent hydrodynamic results. The main 
results of this study can be summarized as follows: 

1. In the proton-rich winds, the proton capture 
can proceed beyond the iron-group nuclei, by by- 
passing the known ^+-waiting point nuclei (e.g., 
^^Ge) via (n, p) reactions, as also suggested by re- 
cent works (Frohlich et al. 2006a, b; Pruet et al. 
2006; Wanajo 2006). These neutrons are con- 
tinuously supplied from the anti-electron neutrino 
capture on abundant free protons. 

2. The neutrino-induced r^^-process leads to the 
production of some light p-nuclei, even with mod- 
erate entropy (~ 50 A^^ A:) that can be found in 
the early winds {tph ~ Is) from the neutron 
star with 1.4 M©. The production of p-nuclei 
is, however, highly dependent on the entropy in 
the wind, which affects the proton-to-seed abun- 
dance ratio Yp/Yh at the beginning of the re- 
process. In the high entropy winds (~ 100 NaIc) 
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with M = 2.0 Mq, the nuclear flow proceeds to 
Z K. 50, resulting in the production of p-nuclei up 
to ^ ~ 130. 

3. The production of p-nuclei is also highly sensi- 
tive to the proton richness of the wind material. In 
the slightly proton-rich winds {Ye ~ 0.50 — 0.52), 
only the lightest p-nuclei (^^Se, ^^Kr, and ^"'Sr) are 
produced. The higher p-nuclei up to A ~ 130 are 
synthesized with interesting amounts as Vg rises 
to ~ 0.56. 

4. The termination of the rp-process is due to the 
increasing Coulomb barrier for the heaviest nuclei 
in the nuclear flow, not due to the exhaustion of 
free protons. The decrease of the temperature be- 
low ~ 2 X 10^ K owing to the expansion of the 
ejecta is also a cause of the termination. The sup- 
ply of neutrons by the neutrino capture on pro- 
tons is suflScient throughout the rp-process phase. 
Hence the dilution of the neutrino flux by the ex- 
pansion of matter is not the principal reason for 
the termination of the rp-process in the current 
wind models. 

5. Few p-nuclei beyond A 130 are produced by 
the neutrino-induced rp-process, owing to the in- 
creasing Coulomb barrier for heavier nuclei. The 
high entropy ('^ 200 A^a^) or highly proton-rich 
(Ye ~ 0.6) wind results in driving the nuclear 
abundances to the /3-stability or even the neutron- 
rich region, when the temperature drops below 
Tg ~ 2. In this case, a variety of neutron-rich 
isotopes, which are generally assign to s-nuclei or 
r-nuclei, are synthesized even in the proton-rich 
compositions. 

6. ^^Mo is produced with interesting amounts in 
the slightly neutron- rich winds (Ye ~ 0.46 — 0.49), 
much more than in the proton-rich winds. This is 
due to the proton capture reactions at T ~ 4 — 
3 X 10^ K, which carry a portion of the nuclear 
abundances from the abundant ^°Zr to ^^Mo. The 
amount of the produced ^^Mo is, however, highly 
dependent on the thermodynamic history of the 
wind. 

7. The neutrino-driven winds can be the origin of 
some light p-nuclei, at least of ^"'Se, ^'^Kr, *^Sr, 
and 92Mo, and hkely of "^Mo, ^^-'^Ru, i^^pj^ 
and -'^'^^'-'^''^Cd, supplemented by the production of 
^^Mo in the slightly neutron-rich compositions. In 
particular, this can be the imique astrophysical 
site responsible for the production of the proton- 



rich Mo and Ru isotopes. The heavier p-nuclei 
{A > 110) may have another origin, most likely the 
0/Ne layers in core-collapse supernovae (Prant- 
zos et al. 1990; Rayet et al. 1995). Note that the 
p-nuclei synthesized in neutrino-driven winds are 
regarded as primary. This is contrast to those 
produced in the 0/Ne layers, which need the s- 
process seed abundances and thus regarded as sec- 
ondary. 

Note that the neutrino-induced rp-process in 
the current study may not contribute to the pro- 
duction of any elements, since the fraction of p- 
process isotopes for a given element is generally an 
order of one percent. Therefore, it would be chal- 
lenging to confirm the occurrence of the neutrino- 
induced rp-process by spectroscopic studies of ex- 
tremely metal-poor stars, which are expected to 
preserve the nucleosynthetic signature from a sin- 
gle (or a few) supernova (e.g., Ishimaru & Wanajo 
1999; Wanajo & Ishimaru 2006b). Analysis of the 
isotope anomalies in some elements found in prim- 
itive meteorites (e.g., Arnould & Goriely 2003) or 
measurements of Galactic cosmic-rays might pro- 
vide us some clues. 

At the end of this investigation, it should be 
cautioned that the nucleosynthesis calculations in 
the current study are based on a schematic for- 
mulation of neutrino-driven winds, which must be 
in fact closely related to the unknown mechanism 
of core-collapse supernovae. In addition, the es- 
timation of the p-nuclei production in each event 
(i.e., the core-collapse supernova that leaves the 
proto-neutron star with 1.4 Mq or 2.OA'/0) was 
based on the Ye distribution during the bubble (not 
wind) phase, obtained from one specific hydrody- 
namic simulation (Buras et al. 2006). Therefore, 
the current results should not be regarded as the 
quantitative predictions that is to be used for the 
Galactic chemical evolution study. 

Nevertheless, the current study provides us 
some notable implications. First, it is likely that 
the neutrino-induced rp-process takes place in all 
core-collapse supernovae to some extent, even with 
a moderate entropy. In this sense, other astro- 
physical sites than standard core-collapse events, 
e.g., coUapsar jets or disk winds formed around 
a black hole, which are associated to gamma- 
ray bursts, can be the astrophysical site for the 
neutrino-induced rp-process. Second, details of 
the nucleosynthesis process are presented for se- 
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Icctcd winds, which wiU be useful for the future 
experimental studies of proton-rich nuclei far from 
the /3-stability. In particular, core-collapse super- 
novae no doubt eject the nuclcosynthctic products 
that contribute to the Galactic chemical evolution 
of heavy nuclei. This is in contrast to X-ray bursts, 
which is unlikely to contribute the Galactic chemi- 
cal evolution. Therefore the experimental estima- 
tions can be easily tested by comparing the nu- 
cleosynthesis calculations to, e.g., the solar com- 
positions. Third, the systematic calculations for 
the wide (but reasonable) ranges of the neutrino 
luminosities and the electron fractions enable us 
to make a meaningful comparison with the solar 
abundances. In particular, a notable aspect is that 
^^Mo has another origin (i.e., neutron-rich winds) 
while other light p-nuclei are mainly synthesized 
by the rp-piocess (i.e., in proton-rich winds) even 
in a single event. This will serve unique con- 
straints to the fluid dynamics of the early super- 
nova ejecta. 
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Fig. 1. — Model parameters {Ru,tph,M) of the 
current neutrino-driven winds arc shown as func- 
tions of I/y. Also denoted are the obtained 
entropies (s) and dynamic timescales (rayn) for 
1.4 M© and 2.0 M© cases. 
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Fig. 2. — Time variations of radius, temperature, 

density, and entropy for selected (odd-numbered) 
wind trajectories with the M = 1.4 model 
{tpb < 2s and tph < 20s for left and right pan- 
els, respectively). The evolution of the neutrino 
sphere is denoted by the grey lines along with the 
staring point for each wind (crosses) in top panels. 
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Fig. 4. — Top: Yg at Tg = 9 (Yei, thin lines) and 
at Tg = 3 {Yef, thick lines) as functions of tp^. 
The solid and dotted lines denote the results for 
M = 2.0 M© with YeO = 0.60 and M = 1.4 Mg 
with Yeo = 0.48, respectively. Bottom: Differences 
of Ye (AYe = Yet - Yef) at tpb = 0.8 s and 4.0s 
as functions of Yeo, for M = 1.4 Mq {dotted lines) 
and 2.0 Mq {solid lines). 
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Fig. 5. Variations of Ye (thick lines) and the 
mass fractions {thin lines) of neutrons pro- 
tons {Xp), a particles (X,^), and heavy nuclei (X^) 
as functions of Tg. The solid and dotted lines de- 
note the results for the wind trajectories W201060 
and W141048 (see the text and footnote 1), respec- 
tively. 
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Fig. 6. — Snapshots of the nucleosynthesis cal- 
culations at t = 0.22 s for the wind trajectory 
of W141060, in which the neutrino-induced reac- 
tions (eqs. [l]-[6]) are turned off {left panel) and 
on (right panel). The temperature decreases to 
Tg « 2 at this time. The abundances are color 
coded in the nuclide chart. The nuclei included 
in the reaction network are denoted by dots, with 
the stable isotopes represented by large dots. The 
abundance curve as a function of mass number is 
shown in the upper left for each panel. 
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Fig. 7. — Same as Figure 6, but for the wind 
trajectories of W201060 {top), W200865 (middle), 
and W200265 {bottom). Left and right panels are 
the snapshots at Tg « 2.0 and ipb ~ 1.0 s, respec- 
tively. Neutrino-induced reactions (eqs. [l]-[6]) are 
included for all cases. 
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Fig. 8. — Nuclear flows defined by equation (9) 
{arrows) and the abundances {circles) in logarith- 
mic scale for the wind trajectory of W201060 at 
Tg = 2.5, for those greater than 10^^. The flows 
by y5+-decays are denoted by off-centered arrows. 
The nuclei included in the reaction network are 
denoted by dots, with the stable isotopes repre- 
sented by squares. 
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Fig. 9. — Time variations of the mass fractions 
{thin lines) of neutrons X„, protons Xp, a. parti- 
cles Xa, and heavy nuclei Xh for the wind trajec- 
tory W201060. The sohd, dashed, and dotted lines 
denote the results with all neutrino-induced reac- 
tions (eqs. [l]-[6]), neutrino capture on free nucle- 
ons only (eqs. [1] and [2]), and no neutrino- induced 
reactions, respectively. Also shown are the average 
mass number of heavy nuclei Ah and the temper- 
ature {thick lines). 
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Fig. 10. Time variations of the mean lifetimes 
for /3+-decay, (n, 7), (p, 7), and (n,p) reactions, 
and neutrino capture on protons for the wind tra- 
jectories of W201060 (top), W200865 {middle), 
and W200265 {bottom), defined by equations (12)- 
(16). The solid and dotted lines denote the life- 
times for Z > 2 and Z > 48 nuclei, respectively. 
For the latter, the values for only Tg < 3 are 
shown, in which these nuclei are synthesized by 
the rp-process. Also shown are the temperature 
and the distance from the center of the neutron 
star {thick lines). 




Fig. 11. — Same as Figure 8, but for the wind 
trajectory of W141048. 
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Fig. 12. — Same as Figure 9, but for the wind 
trajectory of W141048. The mass fractions of ^°Zr 
and ^^Mo are also plotted. 
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Fig. 13. — Mass fractions Xcj of p-nuclei with i^.- 
spect to their solar values (^0, Anders & Grevesse 
1989) as functions of Y^f for M = IAMq (left) 
and 2.0 M© (right) models. Top panels show 
the mass-averaged results (see the text), where 
Yef is taken to be that at tpb — 4 s {L^, = 
2 X 10^^ ergs s~^). The mass fractions of ^^Zn 
and ^°Zr are also plotted for comparison pur- 
poses. Middle and bottom panels show the re- 
sults for the winds with L^, — 8 x 10^^ ergs s~^ 
and 2 x 10^^ ergs s~^, respectively. The histogram 
is the asymptotic Ye distribution p(Ye) of the 
neutrino-processed ejecta during the first 468 ms 
after core bounce for a 15 Mq progenitor star 
taken from Buras et al. (2006). 
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Fig. 14. — Mass fractions of p-nuclei with respect 
to their solar values (Anders & Grevesse 1989) 
as functions of L^, for the selected cases {M = 
1.4 Mq; left, M = 2.0 Mq; right, Yeo = 0.48; top, 
Yeo 0.60; bottom). 
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Fig. 15. — Mass-yg-averaged abundances with re- 
spect to their solar values (Anders & Grevesse 
1989) for M = 1.4 Mq (top) and 2.0 (bottom) 
models as functions of mass number (see the text). 
The abundances smaller than Xcj/Xq < 100 is 
omitted here. The isotopes (after decay) are de- 
noted by open circles (even-Z) and triangles (odd- 
Z). The p-nuclei are denoted with filled symbols. 
The solid lines connect isotopes of a given element. 
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Fig. 16.— Same as Figure 15, but for M = 1.6 M© 
(top) and 1.8 Mq (bottom) models. 
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Fig. 17. — Same as Figure 15, but for L,^ > 
2 X 10^^ ergs (tpb < 4s). Top and bottom 
panels show the results with and without neutrino- 
induced reactions (eqs. [l]-[6]), respectively. 
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Fig. 18. — Same as Figure 15, but for Lj, > 2 x 
10^^ ergs {tpb < 4 s). Top and bottom panels 
show the results for the neutron-rich {¥(./ < 0.5) 
and proton-rich (Ve/ > 0.5) winds, respectively. 
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